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Abstract 

We study the problem of quantum quench across a critical point in a strongly 
coupled field theory using AdS/CFT techniques. The model involves a probe neutral 
scalar field with mass-squared in the range —9/4 < < —3/2 in a AdS^ 
charged black brane background. For a given brane background there is a critical 
mass-squared, such that for < the scalar field condenses. The theory 
is critical when m'^ = and the source for the dual operator vanishes. At the 
critical point, the radial operator for the bulk linearized problem has a zero mode. 
We study the dynamics of the order parameter with a time dependent source J(t), 
or a null-time dependent bulk mass m{u) across the critical point. We show that in 
the critical region the dynamics for an initially slow variation is dominated by the 
zero mode : this leads to an effective description in terms of a Landau-Ginsburg type 
dynamics with a linear time derivative. Starting with an adiabatic initial condition 
in the ordered phase, we find that the order parameter drops to zero at a time 
which is later than the time when (m^ — m^) or J hits zero. In the critical region, 
ti,, and the departure of the order parameter from its adiabatic value, scale with the 
rate of change, with exponents determined by static critical behavior. Numerical 
results for the order parameter are consistent with these expectations. 
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1 Introduction and summary 



An important class of problems in quantum systems is that of quantum quench, where a 
parameter in the hamiltonian varies with time, typically attaining constant values at early 
and late times |^ Starting with some initial state, the problem is to determine the nature 
of the final state. This problem has recently attracted a lot of attention in several areas of 
many-body physics, particularly because of progress in cold atom experiments |l], |2], [s]. 
Among other things, this problem is interesting for two different reasons. The first relates 
to the question of thermalization. Does the system evolve into some kind of steady state 
? If so, is the state "thermal" in any sense ? The second question deals with the situation 
where the quench takes place across a value of the parameter where there is an equilibriium 
critical point. In this case, Kibble-Zurek type scaling arguments indicate that there are 
several physical quantities which are universal and determined by the critical exponents of 
the critical point. Furthermore for two dimensional theories which are suddenly quenched 
to a critical point, powerful techniques of boundary conformal field theory have been 
used in [s] to show that ratios of relaxation times of one point functions, as well as 
the length/time scales associated with the behavior of two point functions of different 
operators are given in terms of ratios of their conformal dimensions at the critical point, 
and hence universal. Another related application of this phenomenon is in gravitational 
physics and cosmology. In this context, this is the phenomenon of particle production in 
time dependent backgrounds, which is relevant to many physical problems ranging from 
black hole evaporation to behavior of quantum fluctuations in an expanding universe. 

There are very few theoretical tools available to study such systems when they are 
strongly coupled. In this note we will explore the use of AdS/CFT correspondence [i] 
- [7| to this problem. AdS/CFT techniques have been in fact used in the past to study 
quantum quench, though not at critical points. In AdS/CFT the boundary values of 
bulk fields become coupling constants, so that this problem becomes that of determin- 
ing a time dependent background with specified initial and boundary conditions. As is 
usual, a quantum problem in the boundary theory becomes a classical problem in General 
Relativity. 

One class of problems involve a vacuum initial state (in the bulk this is pure AdS in 
the distant past) and a change of the boundary value of some field (e.g. the dilaton or the 
boundary metric) over some finite time interval, in a regime where supergravity is always 
valid. In the boundary field theory this corresponds to turning on a time dependent source 

•^Sometimes quench is used to denote a sudden change. We wiU use this word to denote changes with 
arbitrary rates, in particular slow rates. 
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for the dual operator. Under suitable conditions, this leads to black hole formation in the 
bulk [8 10 . The correlators at future time would then be thermal with a temperature 
characterized by the Hawking temperature. The time scale after which this happens 
depends on the nature of the correlators, but turns out to be always smaller than what 
one would expect from a conformally invariant system evolving to a thermal state. Thus, 
in this case thermalization of the field theory is signalled by black hole formation. Another 
class of problems involve a similar setup, but a suitable variation of the coupling which 
prevents black hole formation in the supergravity regime. The coupling, however, becomes 
weak at some time and the bulk string frame curvature grows large, leading to a breakdown 
of the supergravity approximation - thus mimicking a space-like singularity 11 . For the 
case of a slow variation of the coupling it turns out that the gauge theory remains well 
defined and may be used to show that a smooth passage through this region of small 
coupling is possible without formation of a large black hole. Related scenarios appear 
in 



12 and 13 



Many interesting phenomena in gauge-gravity duality, in particular phase transitions, 
can be explored consistently in a probe approximation. In this approximation, a certain 
set of bulk fields can be treated separately and independently of bulk gravity. One set 
of examples are probe branes in a background AdS x 5*. This introduces hypermultiplets 
which in general live on a defect in the original field theory. When the number of such 
branes, Nf is much smaller than the number of flux units which produce the background 
geometry, N^, the dynamics of these branes do not backreact on the bulk metric which 
remains AdS x S. The dynamics is then described by a DBI action in a fixed background 
- the fields in this action are then bulk probe fields. If the boundary value of such a probe 
field is time dependent in a fashion similar to above, black hole formation in the bulk is 
invisible in the probe approximation. However, rather remarkably, thermalization of the 
hypermultiplet sector is still visible. This manifests itself by the formation of an apparent 
horizon (which evolve into an event horizon in some situations) in the induced metric on 



the brane worldvolume 14 , 15 . Fluctuations of the brane feel this induced metric and 
their correlators have thermal properties |^ 

None of the above works study quench dynamics near a phase transition. In this 
paper we take the first step towards doing that, in a simple model which displays an 



equilibrium critical point in a probe limit 19 . The model involves a neutral scalar 



field whose mass lies in the range —9/4 < < —3/2 in the background of a charged 
AdS^ black brane. The overall coupling is chosen to be large so that the scalar can be 
treated as a probe. For a given brane background (i.e. for a given temperature and 



*For related phenomena see 
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charge density in the dual field theory), and for a vanishing non-normalizable mode of 
the field (corresponding to a vanishing source of the dual operator) there is always a 
critical mass below which the trivial solution is unstable. In this regime the stable 
solution is nontrivial, signifying a nonzero expectation for the dual operator in the dual 
theory, = with a vanishing source is then a critical point, which happens to have 
standard mean field exponents. This setup is in fact quite similar to models of holographic 



superconductors and fermi surfaces 20 - 26 . However, here the scalar field is neutral so 



that the dynamics of the gauge field is frozen. This model has been argued to model 
phase transitions in antiferromagnets. As is well known, in this mass range there are two 
inequivalent quantizations. Our analysis is restricted to the conventional quantization 
of the bulk scalar, where the mode which vanishes slower at the boundary is treated as 
the source in the boundary theory. A similar analysis can be easily performed for the 
alternative quantization. 

We will study the time evolution of the order parameter starting from adiabatic initial 
conditions in the ordered phase. The passage through the critical point is done in two 
ways - (i) by keeping the bulk mass at its critical value and changing the leading term of 
the near-boundary expansion of the bulk field and (ii) by changing a bulk mass parameter 
with a vanishing leading term. The field theory meaning of (i) is clear - this corresponds 
to turning on a source dual to the bulk field : the critical point appears when the source 
vanishes and the mass is at its critical value. The meaning of (ii) is less direct - changing 
a bulk mass means changing the anomalous dimension of the corresponding operator in 
the field theory. As discussed in 19 , one way to do this is in fact coupling this scalar 
field to another field and changing the source for this other field. 

Away from the critical point, adiabaticity holds for a sufficiently low rate of change. 
We will show that the adiabatic corrections are proportional to linear time derivatives, 
even though the equations of motion involve only second order time derivatives. This is a 
well-known phenomenon and happens because of dissipation implied by the presence of a 
horizon. In the bulk, we need to impose regularity conditions at the horizon, which is easily 
done using ingoing Eddington-Finkelstein coordinate |27||28] u as time, u coincides with 
the usual time coordinate t on the boundary and d/du\r = d/dt\r, though d/dr\u 7^ d/dr\t. 
The bulk equations now involve first order derivatives S/Sm]^ As expected, adiabaticity 
fails as one approaches the critical point at a time which we choose to be t = 0. The 
appearance of a first order time derivative is related to the fact that the dynamical critical 



29 



''Note that the analysis can be done using the usual time coordinate as well, as originally done in 
but one has to take special care of the behavior of the solution near the horizon. The result is of course 
the same. 
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exponent of similar models turn out to be 2; = 2 30 . 

When adiabaticity fails, the system enters a scaling regime. A useful way to understand 
the dynamics is to decompose the bulk field in terms of eigenstates of the radial operator 
which appears in the linearized problem around the equilibrium solution. For any nonzero 
temperature of the background black brane, this operator has a zero mode which is regular 
at the horizon exactly at the critical point. We will show that in the critical region, the 
dynamics for small v is dominated by that of this zero mode. The other modes remain 
adiabatic. The scaling properties can be then understood in terms of a Landau-Ginsburg 
type dynamics of the zero mode with a linear time derivative. The order parameter 
remains nonzero even when the instantaneous value of the coupling reaches the equilibrium 
critical value and first drops to zero at a time ~ f ~" where v denotes the rate of change 
of the coupling in the critical region. The departure of the value of the order parameter 
from the instantaneous equilibrium value also scales as . The exponents a and /3 are 
essentially determined by the static critical exponents. 

Beyond the critical region, however, the other modes become important. In fact, at 
very late times the bulk equation may be approximated by its linearized form, so that the 
solutions are nothing but the usual quasinormal modes around the background black hole. 
As usual, the imaginary part of the lowest quasinormal mode frequency then determine 
the decay of late time behavior. 

We also present preliminary results of a numerical solution of the time dependent 
equations for the bulk field and extract the time dependence of the order parameter. The 
results are consistent with the above picture. Our numerical results are not yet accurate 
enough to verify the scaling behavior. However this should be possible with further work. 

In Section 2 we discuss spatially homogeneous classical quench in a Landau-Ginsburg 
type model as a prelude. In Section 3 we review the main faetures of the equilibrium phase 
transition. In Section 4 we discuss the failure of adiabaticity in quantum quench across 
this transition. We then show that in the critical region the dynamics for an initially 
slow quench is dominated by that of the zero mode, and discuss its relationship of the 
dynamics to the Landau-Ginsburg model of Section 2. In Section 5 we present the results 
of a numerical solution of the bulk equation of motion and verify the expectations in 
the previous section. Section 6 contains some concluding remarks. The three appendices 
provide details of derivation of the some of the results contained in the main text. 
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2 Classical quench in a 0+1 dim model 



As a prelude to the main discussion, in this section we will study classical spatially 
homogeneous quench in a Landau-Ginsburg model, which is described by the equation, 

'^"'^+m2(t)0 + 02/3-1 + j(^) = (2.1) 



As discussed above, holography converts the problem of a quantum quench to a problem of 



classical quench in one more dimension - which is the motivation behind studying (2.1). 
In a later section we will show that this simple model describes some of the essential 
features of the dynamics in the critical region of the holographic phase transition which is 
the main subject of this paper. Note that we do not have any noise term (which is what 
is needed for many other applications). 

When and J are time independent, this model has a critical point at = J = 0. 
Here we have kept generic a and (3 for convenience of later discussion. In most conservative 
models, i.e. where energy is conserved, a = 2. However as we will see a = 1 is more 
appropriate for the holographic transition which we will discuss. In the following we 
will therefore use this value. We will also concentrate on the /3 = 2 case, i.e. a quartic 
interaction. 



2.1 J Quench 

Let us first discuss the case where J{t) is time varying, but a which is independent of 
time and positive. 

For a time independent J we have, 

m'^cj) + (P^ + J = (2.2) 

Denote the solution of this static equation by (f)o{J,m). We then introduce a J{t) which 
is slowly varying with time. The idea is to start with adiabatic initial conditions at some 
early enough time and study the time evolution of the order parameter. To study the 
adiabatic expansion we write, 

= 0o(J(t),m) + e0i(t) + --- (2.3) 

Here e is an adiabaticity parameter introduced by performing a scaling t — )■ t/e so that 
any dt comes with a factor of e. The ellipsis denote higher order terms in e. To lowest 
order in e, the equation governing (pi is, 

0i(m2 + 30^(t)) = 0,(t) = j(t)^ (2.4) 
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The adiabatic expansion is good for a finite m,(j)Q for a sufficiently slowly varying J{t). 
The adiabatic expansion fails when 



~ (1)0 (2.5) 

For a sufficiently slowly varying J{t), this is possible if both and 0o are small, i.e. when 
we are close to the critical point = J = 0. Approaching the critical point along the 
direction = we have 0o(^) ~ {^■Ji^))K and 
of adiabaticity , 



2.6 



leads to the condition for breakdown 



j~j3. (2.7) 
If the quench is linear near the critical point, that is J(t) ~ Jovt for small J, we get 

vtl ~ 1 (2.8) 

as the condition for breakdown of adiabaticity. 

When adiabaticity breaks down, the system enters a scaling region. In the scaling 
region we use Jit) = J^vt + Oit^). We rescale the field as: (p — )■ a0 and t — )■ (31. Choosing 
0^(3 = 1 and f ^ = 1, we get a v independent equation for The quadratic or higher 
order temporal contribution in J{t) only gives a sub-leading contribution in Jqv in the 
scaling limit. Hence we have a scaling solution. 



1 ~ 



(f){t) =vy^(f){v5t). (2.9) 

This scaling solution leads to an estimate of the magnitude of fluctuations {6(f)(0)) ,i.e. the 
departure of (j) from the equilibrium value at t = 0. Outside the critical region 50(0) ~ 
whereas in the critical region 50(0) ~ v^. This also defines the "zero crossing time"(t^) 

2 

defined by (f)(t^,) = 0. We have ~ f ~5 at the critical point. Which diverges in the f — 
limit. Out of criticality actually approaches a constant in the f — )■ limit. 

It is instructive to study the dynamics of the order parameter for a profile of J{t) 
which becomes a constant at early and late times and behaves linearly with time in the 
critical region, for example 

J(t) = Jo tanh(f t) (2.10) 
with = for all times. The numerical solution of the equation of motion with an 



adiabatic initial condition at early times is shown in Figure (2.1) for a typical value of v. 

Numerically we see that at sufficiently early times the order parameter tracks the 
adiabatic solution and overshoots the adaiabatic solution as one approaches the critical 
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Figure 1: The order parameter as a function of time for a J{t) = tanh(ft) at = 
with V = lO^*^'^, 10^^, 10^^'^, 10^^ (from the bottom on the left). The adiabatic solution 
(dashed) is also shown as a comparison. 



region at t = 0, falling down to zero at some positive time t^. At later times, the solution 
again approaches the adiabatic solution. The numerical solutions for different values of 

1 2 

V are consistent with the behavior 50(0) ~ t>5 and tj, ~ f "s as predicted by the scaling 
analysis. 

For profiles of J(t) which behave as some nontrivial power of t in the critical region, 
i.e. J{t) ~ (ft)" a similar analysis leads to 



n ~ An 
(f){t, V) = V 2"+3 (f)[tV 2"+3 , 1) 



(2.11) 



leading to the scaling properties 



+3 



__2n_ 



(2.12) 



If there are more than one interaction terms, e.g. 0^^"^ present in the equation of 
motion, then the lowest order interaction (/3o) dominates the critical dynamics. In this 
case we have, 

n ~ (2/3o-2)n 

(pit.v) = t,(2/3o-2)(n+l) + 10(tt,("+l)(2/3o-2) + l). (2.13) 



2.2 Mass quench 

In this section we will study classical quench across the critical point along the line J = 
by changing from negative to positive values. For static < 0, the model has a 
non-trivial time-independent solution 0o = V—rn^- The static model has a second order 
phase transition near = 0. Now we introduce a time dependent mass function which 
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asymptotes to some constant value of < at early times and then rise to a constant 
value > at late times, e.g. 



m'^{t) = -mltanh{vt) (2.14) 

Once again we start with adiabatic initial conditions at some early enough time and study 
the time evolution of the order parameter. To study the adiabatic expansion we write, 

= 0oW + e0iW + --- (2.15) 

where, 

0o(t) = V-m2(t) (2.16) 
To lowest order in the adiabatic expansion 

^ (2.17) 



2m2 



Using the fact that 0o ~ y—iTT-'^), we get that the adiabatic expansion breaks down as 
01 ~ 005 i-G- when 



1. 2.18 



As expected this breaks down as we approach the equilibrium critical point. If (— m)^ ^ vt 
near the critical point breakdown of adiabaticity happens when, 

(2.19) 

When adiabaticity breaks down the system enters a scaling regime. In this regime, we 
can approximate m?{t) = —mfvt. Here the scaling relations are 

(j){t,v) ^v-^4){tv^,l) (2.20) 

However, if the solution does not cross a phase transition point then fluctuation around 
a static solution is suppressed exponentially in ^. It is also to be noted that if m{t) is 
sufficiently slowly varying, the exact functional form of m(t) is not needed. The only 



information we use to derive (2.20) is that m{t) is a linear function of time near the 



phase transition and m{t) — )■ m2 as t — ?■ oo. 
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3 The equilibrium phase transition in the holographic 
model 



The model of 19 has a neutral scalar field r, x) in the background of a charged AdS^ 
black brane. The lagrangian is given by 

The background metric is given by (in RAds = 1 units) 

ds^ = [-r2/(r)rft' + rV] + (3.22) 



where 

The associated Hawking temperature is then given by 



/(r) = [1 + ^ - i±|^] 0<r^<l (3.23) 



T = -^{1-V) (3.24) 
Airro 

In the following we will replace r — > rrg. 

One should add the Einstein- Hilbert action and a cosmological constant term to (3.21 ). 
In the limit of large A, however, the field can be regarded as a probe field since its back 
reaction to the metric via Einstein's equations is small. This is the approximation we 
adopt. 



In 19 it was shown that when the mass lies in the range 

-1<^'<-1 (3-25) 

there is a critical phase transition at some value of T = Tc{m) when the source to the 
dual operator vanishes. Conversely, for a given T there is a value of = where the 
theory is critical. 

The upper limit in (3.25) is the BP bound for the near-horizon AdS2 geometry which 
appears in the extremal [rj = 0) metric. (Note that the AdS scale for this infrared AdS2 
is given by l/V^ in our units). The lower bound is the BP bound for the asymptotic 
AdS^. Pield configurations which are translationally invariant in the x directions satisfy 
the equations of motion 

^^-J^^t + ^r(rV(r)9,)]0 - - 0' = (3.26) 
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Near the AdS^ boundary the asymptotic behavior of the solution to the hnearized equation 
is of the form 

0(r) = J{t)r-^'[1 + 0{l/r^)]+ <0> (t)r-^+[l + 0{l/r^)] (3.27) 
where A is given by 

A± = ^ ± ^^' + 1 (3-28) 

In the range of masses of interest, both the solutions are normalizable, so that there is 
a choice of quantization. The standard quantization considers the coefficient J{t) is the 
source in the dual field theory and < O > (t) then gives the expectation value of the dual 
operator. In the alternative quantization the expectation and source change the role. 

Consider first the linearized problem, ignoring the cubic term. By a standard change 
of coordinates to tortoise coordinates p and a field refinition to 

* = -^ #r.<) = ^ (3.29) 

The horizon is then at p = oo and the boundary is at p = 0. the equation becomes ( [24j) 

- dh = -dlx + Vo{p)x = V,x (3.30) 

with 

Voip) = f{r)[{m' + 2) - ^ + (3.31) 



where in Vo(p) we need to express r in terms of p using (3.29). 

For solutions of the type x ~ e~*'^*, equation (3.31) is a Schrodinger problem in a 



potential V{p). The potential goes to zero at the horizon p = — oo and behaves as ^"^^^^-^ 
near the boundary p = 0. Note that for a brane background at any finite temperature, 
/(r) ~ (r — 1) near the horizon, while p ~ — log(r — 1) so that Vq ~ as we approach 
the horizon. In contrast, for the extremal background /(r) ~ (r — 1)^ while p ~ l/(r — 1) 
so that Vq ~ 1/p^- This makes the analysis for the extremal background rather subtle. 
In this paper we will work with the non-extremal case. 

Due to the vanishing of the potential near the horizon the time-independent problem 
has a continuum of modes starting at w = 0. The behavior of the solution near the horizon 
p = —oo is of the form e~*'^'^*=''''\ the two signs representing ingoing and outgoing waves. 
We are eventually interested in solving the problem with ingoing boundary conditions at 
the horizon. 



^When we turn on J(t), it is a valid concern whether we will be able to neglect the non- linear term 
near the boundary. This can be done as long as A > or to^ < 0. 
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1.0 - 



Figure 2: Profile of (rescaled) zero mode (x ^~(j){x), wliere x ~ ^)) witli —2.235 
and 7] = 0.98. 

Generically for any value of there are quasi-normal modes with complex frequencies. 
These are normalizable modes with ingoing boundary condition. Unlike the continuum 
discussed above these modes form a set of disconnected poles in the complex frequency 
plane. However, normalizable bound states appear when is sufficiently negative : this 



is what has been shown in 19 



There is a critical value of = when a zero energy bound state appears, which 
vanishes in an appropriate fashion at the boundary and is in addition purely ingoing at 
the horizon. In the complex frequency plane some quasinormal mode(s) hit the origin at 

= m^. This critical mass is = — | when r] = 1 and decreases with decreasing i] or 
increasing temperature. The existence of such a bound state for = —2 has been shown 



by variational methods in 21 . A similar demosntration should be possible for other values 
of m^. We have numerically verfied the existence of this zero mode. A typical plot the 
zero mode is shown in in Fig [2j Note that the leading large-r behavior has been factored 
out to ensure that the source term indeed vanishes at the boundary. 

The presence of a bound state in the Schrodinger problem means that the solution 
X = of the full nonlinear problem is unstable, and one has to look for other nontrivial 
solutions. In [19| it was shown that such a stable nontrivial static solution (f)o{r) exists in 



the range of masses given by (3.25), both for the standard and alternative quantizations. 
This means that the expectation value of the dual operator is nonzero, i.e. the field 
condenses, and = m^(T) is a critical point . One can approach this critical point 
either by fixing the mass and tuning the temperature, or by fixing the temperature and 
tuning the mass. The latter may appear to be somewhat strange from the point of view 
of the boundary field theory. However, it is only the IR mass which is relevant - and one 
may imagine obtaining this by coupling to another field \l/ by The expectation 
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value of \I/ which may results from changing a coupling in the boundary field theory then 
renders a mass to 



19 



The critical behavior has standard mean field exponents at any finite T. If the operator 
dual to the field (p is O and the source is J then an analysis identical to that presented 
in 



19 leads to (for — — )■ 0+) 

d<0> 



<o> 



(m 



■m 



2U/2 



dJ 



I j=o 



J^/^ (3.32) 



Exactly at zero temperature the phase transition is of BKT type and the order parameter 
depends exponentially 



< O >j=o~ exp 



2^r 



(3.33) 



The zero mode of the linearized operator will play a key role in the following. This zero 
mode is known to exist for any non-zero temperature. At zero temperature the situation 
is less clear. For a similar zero temperature system involving the D3-D5 system 32 such 



a zero mode does not exist 33 



4 Quench across the critical point the breakdown of 
adiabaticity. 

Our aim is to study quench dynamics across this equlibrium phase transition. In a quench 
situation, some parameter of the boundary theory is made time dependent. The first issue 
we will study is the breakdown of adiabaticity as we approach the critical point. First, 
we argue that away from the critical point, the order parameter tracks the changing 
coupling adiabatically. We will then determine the regime near the critical point where 
adiabataicity breaks and show that in this region the order parameter and the time scale 
are proportional to some power of v pretty much like the model studied in section 2. 



4.1 Quenching the source 

As discussed above, a simple way to quench across the critical point is to remain at the 
critical value of the mass and consider a time dependent source for the dual operator in 
the boundary field theory. This corresponds to a nontrivial boundary value for the bulk 
field. Our discussion is in the conventional quantization of the theory, i.e J{t) in (3.27) 
is treated as a source. The case of alternative quantization should follow along identical 
lines. 
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It is well known that to study low frequency modes in the background of a black brane 
it is convenient to use ingoing Edddington-Finkelstein coordinates, 



u = p — t, 



P 



(4.34) 



where p is defined in (3.29). In terms of these coordinates the equation of motion 3.26) 
becomes 



where 



-2dud,x = -d'^X + V{p,x). 
V{p,x) = Vo{p)x+^x'. 



(4.36) 



This equation has to be solved with the boundary condition that the field is regular 
at the horizon, which at the linearized level is equivalent to requiring that the waves are 



purely ingoing at the horizon 27,28 . At the linear level the analysis can be of course 



performed in (t, p) coordinates as well. However, in this case one needs special care to 
extract a leading singular piece before performing the low energy expansion |29|. At the 
full nonlinear level this procedure is possibly rather involved, but imposing regularity in 
EF coordinates remain simple. 



We need to solve (4.35) with the condition 



x{u,p)-^p ^^'^'J{u) as p — > 



(4.37) 



where A± are defined in (3.28). We have considered the source to be a function of u 



rather than a function of t. On the boundary p = p and t are the same (and quite 
generally d/du\p = d/dt\p), so that this does correspond to a time dependent source in 



the boundary theory. However, as we will soon see, it is convenient to use (4.37). This 
form makes bulk causality explicit. 

To perform the adiabatic expansion, let us decompose the field xip^u) as 

Xip,u) = xiip,u) + Xsip,u) (4.38) 

Where XiiP^u) = J{u)p^^^^' and Xs{p,u) ~ p^^^^+ as p — )■ 0. We get, 

[-dl + Vo{p)]xs + ^ix's + ^Xixl + ^xlXs) = -[-dl + V,{p)]xi - (4-39) 

-2dudpXi - 2dudpXs 

If xi{p,u) = Xi{p) is time independent we have a static solution Xs{p,u) = Xo{p)- This 
has been discussed in the previous section. When we turn on a time dependent J{u) 
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which is slowly varying, one expects that the time dependent solution of (3.26) may be 
constructed in an adiabatic expansion. To understand that we write, 



Xs{p,u) = xoip, J{u)) + e Xi{p,u) + •. 



(4.40) 



Here e is an adiabaticity parameter which keeps track of the adiabatic expansion. If we 



scale M — > w/e, each u derivative is of order 0(e). The ellipsis in (4.62) are higher order 



terms in e. The idea then is to insert (4.62) into the equations of motion and obtain 
equations for xi, X2, • " " order by order in e. To the lowest order one gets 



{[-dp + Vo{p)] + ^(3xo + 6x^X0 + 3xf)}xi = +2dud,xi + 2dud,xo (4.41) 



In deriving (4.41) we have ignored any time dependence of xi? as is appropriate in an 



adiabatic expansion. Generically, at each order in perturbation we may cast the adiabatic 
problem as, 

V<^'^'>Xn = Jn{n,p), (4.42) 

where Jn{u) is a source term contains time derivatives and determined by the lower order 
equations in an adiabatic expansion. 

As discussed before, near the black hole horizon the differential operator becomes a 
plane wave operator and the potential part vanishes. Since the potential Vo(p) goes to 



zero at the horizon p = oo and the terms which involve xo, Xi (4.41 ) also vanish at the 



horizon (because of the overall factor of /(r)), the operator Vp has a continuous spectrum 
which begins at zero. It is, therefore, a non-trivial fact that there is adiabaticity away 
from the critical point. We will now show why this is so and how adiabaticity breaks 
down near the critical point. 

These equations are most conveniently solved by Green's function method. 



Xn = / dp'G^"\p,p')Mu,p'). 



(4.43) 



where G^"'\p,p') is the Green's function of the operator T>p^^ with the boundary conditions 
G(0,p) = and G{oo,p) is regular. 

Let us concentrate on the lowest order terms in the adiabatic expansion, n = 1. All 
higher terms may be found out by similar procedure. The Green's function is given by, 

1 



G^'\p,p') 



" Wixi,X2] 

1 

W{xi,X2) 



■Xi{p')X2{p), P<P 
X2{p)xi{p), P>p', 



(4.44) 
(4.45) 
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where xi ^^id X2 are solutions of homogeneous part of eqn (4.42) satisfying appropriate 
boundary condition at the horizon p = oo and the boundary p = respectively, and 
W{xiiX2) is the Wronskian which is independent of p in this case. We have normalized 
Xi(p) and X2(p) in such a fashion that xi = 1 at the horizon and X2 p^^^^~ near the 
boundary. 

Given that we are looking for a solution which is regular at the horizon (p — )■ oo), xi 
is a constant at the horizon and so is the Green's function. Nevertheless, the integral in 



eqn (4.43) is is finite. This is because regularity of the functions xi and Xo mean that 
drXhdrXo are finite at the horizon. Since (r — 1) ~ e"'', this implies that J{u,p) ~ 
dpixi + Xo) ~ 6xp(— p), damping out any divergence in (4.43) from the large p region. 



In general, near the horizon X2 can be expressed as a linear combination of a regular 
and irregular solution, i.e X2(p oo) = ap + b. It can be easily checked that W{xi^ X2) = 
a. Hence, 

Xi{p,u) = / dp' [e{p - p')xi{p)X2{p) + {p^ p')] {2d^d,Xo + 2d^dpXi) (4.46) 



This is finite as long as a is finite, and small for sufficiently small duXo and duXi- 

At the phase transition point, however, the operator V^p^ has a zero mode which is 
regular at the horizon. This is because at J = the factors involving xo and xi in the 



left hand side of (4.41) vanish, and the operator is identical to the operator acting on 



the linearized small fluctuations at = around the trivial solution xo = 0, i.e. the 



operator Vp which appears on the right hand side of (3.30). We know that this operator 
has a zero mode which is regular at the horizon and vanishes as p~^~*"^~ at the boundary 
p = 0. This means that at this point a = 0. Therefore, the first adiabatic correction 
diverges. 

To derive the precise condition for adiabaticity, it is sufficient to look at the operator 
V^p^ for small J{u). For small J{u), the leading departure from the critical operator comes 
from the term which is proportional to Xo ~ J'^^^- Thus we can use perturbation theory in 
J to estimate a oc J^/^. As argued before, Xo ~ while xi ~ -J- Hence the leading 

divergence in xi can be estimated as 

Xiip,u)r~.:l^ = J-^/^j. (4.47) 

The breakdown of adiabaticity happens when, 

Xi(p,m)~Xo (4.48) 
^ j ~ J^/^ (4.49) 
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exactly as in the + 1 dimensional model (equation (2.7). Assuming a linear quench, i.e 

J{u) ~ vu, we get 



5 



vu'2 ~ 1 (4.50) 
as the condition for breakdown of adiabaticity. 



4.2 Scaling near the phase transition 



In this section we argue that the critical region dynamics is dominated by that of the 
zero mode for (initially) small quench rate. We will demonstrate this for the case where 
J{u) = vu in the critical region. Extension to non-linear quenches is straightforward. 

1 ^ 2 _ 

To understand this let us rescale x ~^ v^Xs^u — )■ v~^'U.. Assuming a linear quench, 
i.e. xi ~ 'vuxiip), we have at the leading order in v, 



[-dl + Vo{p)]Xs + vl l^iXs? + u[-dl + Vo{p)]xi + 2d^d,Xs] + 







(4.51) 



The ellipsis denote terms which contains higher powers of v. 

Now let us expand the sub-leading part of the scalar field in eigenfunctions of the 



operator Vp (defined in equation (3.30) at the critical point 



Xs{p,u) 



ak{u)xk{p)dk 



where Xk satisfy 



(4.52) 



(4.53) 



'P'pXk = l-dl + V,%p)]xk = k'xk 
where denotes the potential in (3.31) at = m^. the eigenfunctions Xk{p) are delta 



function normalized. The eigenfunctions obey the condition 

Ump^oip'^'Xkip)] = 



(4.54) 



In terms of the eigen-coefiicients ak{u) the equation (4.51) becomes, 
fc^afe + v^' ( uJk — I bkk'dudk'dk' — I ak'dk"dk"'Ck,k' ,k" ,k"'dk' dk" dk" 



where, 



Xkip)[-dl + Vo{p)]xi dp 



bkk' — J dpXkdpX'k ! Ck,k',k",k"' — J 



dpXkXk'Xk"Xk"' 



fir) 



(4.55) 
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The equation (4.55) suggests that there is a solution in a perturbation expansion of 

~ 2 

dk{u) = S{k)^o{u) + V5fik{u) H , 



powers of f 5 . Writing 



we get an equation for 



uJo- boo^Uu) - Cooooeo(w)' = (4.56) 



and a set of equations for 7]^. 

fikiu) = (^uJk - bko^loiu) - Ckomloiuf^ (4.57) 

= {u{Jk - Jo) - (bko - boo)-^Uu) - (Cfeooo - Coooo)fo(w)') ) (4.58) 

In the last line of the above equation we have subtracted the equation for Xo(w) to make 
the — )■ limit explicit. 

If the values of k were discrete, the k = mode would clearly dominate the dynamics 
for small v. However, is a continuous parameter starting from zero - this naively indicates 

2 

that the perturbation expansion in f s could break down. Whether this happens or not 



depends on the behavior of the numerator on the right hand side of (4.58). 

It turns out, however, that all the terms like {Jk — Jo) ~ fc^ as A; — )■ 0, rendering the 
perturbation expandion in ft well defined. This is shown in Appendix (|c|. 

Thus, in the critical region the solution Xs is dominated by the zero mode piece 
proportional to Xo- terms of the original variables, we have finally 



Xsip, u) ^ v^^{v^>u)xoip) + y Vkiv^>u)xkip)dk (4.59) 
Since the equation for is identical to the equation for the quantity in the toy model 



(2.1), we are led to scaling relations identical to the + 1 dimensional model. The one 



point function of the dual operator is given by 

< O > (n) ~ Lim,^o[p'"^+Xs(p, u)] (4.60) 

Therefore, the time dependence of O is governed by the same scaling behavior. 

The above analysis can be carried out for away from the critical value as well. In 
this case, the terms like JT^ ~ as — > 0, so that the expansion in f is not valid. This 
means that away from the critical point, all the modes are important. 
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4.3 Mass quench 

Another way to quench across the critical point is to consider a time dependent with no 
source, keeping the background geometry fixed. The mass function is taken to asymptote 
to some constant value of < at early times and then rise to a constant value 
> ml at late times, e.g. 

m'^{u) = ml + {ml — m\) tanh(wM) (4.61) 

Note that we have a mass which depends on the EF coordinate u rather than the time 
coordinate t. The motivation for this is as follows. As mentioned earlier, we may consider 
the bulk mass term as arising from a coupling of the scalar field with some other field. 
This latter field acquires an expectation value in response to a source which leads to a 
mass of our field (f). In this scenario, a time dependent source of the second field would 
lead, by causahty, to an effective mass-squared which has the form m'^{u)G{p) + ■ ■ ■ , 
where the ellipsis denote terms which involve u derivatives of m{u). The function G{p) 
vanishes at the boundary and is regular at the horizon. To leading order we therefore get 



an equation like (4.35) with the potential Vq modified by replacing — )■ m'^{u)G{p). 

The analysis in the previous sections can be now easily repeated for this problem. 
Below we outline the main steps and results. 

Let us consider the case of a time dependent mass which asymptotes to a sufficiently 
negative value, mf. If the mass remained constant at there is a nonzero static solution 



Xo('", ""^i)- Now once we make the the mass time dependent we solve (4.35) adiabatically 



in du- For a mass which is slowly varying one expects that the time dependent solution 



of (3.26) may be constructed around an adiabatic solution, 

X(p, u; m{u)) = Xo{p, m{u)) + Xi{p, u) + (4.62) 

Xo(p; ^) is the static solution obtained by ignoring the time dependence of the mass com- 
pletely. Here e is an adiabaticity parameter which keeps track of the adiabatic expansion. 

Following similar arguments like J quench case , the divergence in xi can be estimated 
to yield 

Xi{p,u) ~ ^duiV^^'i^)) (4.63) 

where 5m'^ = [ml — m^). The condition for break down of adiabaticity for linear quench 
{5m^{u) ~ vu near the critical point) is similar to + 1 dim case, 

tv^ ~ 1 (4.64) 
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When adiabiticity breaks down, the system enters a scahng regime just hke the + 1 
dim case. Following arguments similar to what discussed in the previous section we find, 

X(p, u) ^ v^q{v^u)xo{p) + J bk{v^u)xk{p)dk (4.65) 

The first term defines the scaling relation. The expectation value of the scalar field in the 
scaling regime has the following properties, 

< 0{t) >=< 0{tv'2) > v'^ (4.66) 



5 Numerics 



In this section we will solve the equations (3.30) numerically, although our analysis is quite 



preliminary. Mainly we would like to understand how the system behaves away from the 
regime of slowly varying quench. It is generally hard to solve PDE's numerically. However 
we will use the convenient method of lines. In this method the PDE is at first discretized 
only in the r direction. This gives rise to a set of ODEs. This set of ODEs is then treated 
in a ODE solver. We fix the boundary condition at the two ends, at the boundary and 
at the horizon. We use a regularity boundary condition at the horizon and appropriate 
standard quantization condition at infinity. The initial radial profile of the scalar field 
and its derivative are given at t = 0. In most cases the initial value scalar profile would 
be a solution with = m^(to)- Hence < m^. Then we will change to some other 
value higher than the critical mass. The first point to note is that at late times , i.e. 
out of the critical region, any finite energy excitation is gradually sucked up by the black 
hole. In fact at late times, the equation can be approximately linearized and the solution 
in a black hole background can be expressed as a sum of quasi normal modes. Hence, the 
quasi normal frequency with smallest imaginary part determines the late time decay of 



the scalar field in a black hole geometry 36 . We find a similar decaying behavior even 



in the full non-linear theory (Fig [3]). Note that the scalar field falls in the black hole as 
expected. 

A time dependent m{t) changes the scalar field profile near the boundary. As the 
potential vanishes near the black hole horizon, changing the value of m{t) does not have 
an immediate effect. Hence a time dependent m{t) creates disturbances near the boundary. 
The disturbances then propagate towards horizon along a light cone. 
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(a) The plot is for (f>{x,t), {x ^ ^). It is (b) Plot of boundary expectation 

to be noted how the disturbance propa- value < O > (t). 

gates following a light cone. 

Figure 3: A typical response under a continuous but fast change of m^(t) ^ —2.1825(0.9 — 
0.1tanh(3t)). Here t] ^ 0.991 and ^ -2.169 at t = 0. 

6 Outlook 

In this work we have taken a first step in using AdS/CFT methods to study quantum 
quench across holographic critical points. The particular critical point we study is a 
mean field transition at any finite temperature, and it is not surprising that the critical 
behavior can be obtained by a Landau-Ginsburg type model of an appropriate mode. 
It may appear surprising that the dynamics involves a first order time derivative, even 
though the underlying theory is Lorentz invariant. As discussed above, this is a feature of 
non-zero temperature which manifests itself as a dissipative horizon in the bulk. However, 
similar techniques can be used to study other interesting non-mean field theory transitions. 
In particular, the zero temperature limit of the transition which we study is of Berezinski- 
Kosterlitz-Thouless type. Other examples of such non-trivial transitions are discussed 
in 32 . The situation is more subtle in this case : it would be worthwhile to study quench 



in these systems. It would be also interesting to study dynamics of holographic transitions 



associated with black branes whose IR geometry is of Lifshitz type 38 . 

As discussed above, AdS/CFT maps quantum quench in a strongly coupled theory to 
a classical quench in one higher dimension. The dynamics is completely deterministic, as 
opposed to e.g. Langevin dynamics - the effect of a non-zero temperature being accounted 
by the presence of a black brane in the bulk. Such systems have also been studied in the 
past [31 1 . In these studies, however the interesting questions are rather different. They 
involve, for example, long time behavior of fluctuations and correlation functions when 
averaged over initial conditions. These questions can be addressed in our framework as 
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well - they would correspond to corrections of the dual field theory and are left to 
future work. A related direction is the study of entanglement entropy during a quench 



which has been studied holographically in other contexts in 9 - 10 



Finally, it would be interesting to extend our considerations to inhomogeneous quench. 
In this case one interesting question relates to defect formation as one crosses the phase 
transition. Some aspects of this question has been studied numerically in 37 . It would 
be interesting to see if expected scaling properties continue to hold, and if there is any 
analytical insight in terms of the eigenmodes of the radial operator. 



A a = 2 quench 

In this section we will study quantum quench in a simple relativistic Landau-Ginsburg 
model. The relevant lagranian is 

In our language this is a system with a = 2. Quench in a holographic model without any 
black hole horizon, e.g. a quench in global AdS, is expected to obey eqn (22-1) as an 
effective description. This is suggested by the fact that in this latter situation the lowest 
order adiabatic correction is proportional to the second derivative of the source, as is 



consistent with lack of dissipation (see e.g. the last reference in 11 . This kind of model 



also appears in dynamics of phase transitions in cosmology 39 A thorough discussion on 
those issues is left to the future. 



The lagrangian (A.67) leads to the equation of motion 

+m^{t)(l) + ^^ + J{t) = (A.68) 

When and J are time independent, this model has a critical point at = J = 0. We 
study quench by seting = and considering a time dependent J{t) = tanh(ft). For 
small V, breakdown of adiabaticity occurs when, 

tv'i ~ 1 (A.69) 

The scaling solution in the critical regime is given by, 

= (A.70) 



^We have rescaled fields to set the self coupling to 1. 
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As discussed in the main text this scahng relation fixes 6(f)(0) and zero crossing time t*. 

If we continue the quench we may ask how much excess energy (or the equivalently 
the amphtude of oscillation) remains in the system at late times. If the adiabatic approx- 
imation were valid, the excess energy would be suppressed exponentially. However, here 
the system goes through a critical point where adiabaticity is broken and at late times 
the energy turns out to scale as a power of v. 



B k dependence of coefficients : A toy model : 

The purpose of the next two appendices is to justify the perturbation expansion in f 
used to derive the critical region scaling in section 4.2. The idea is to examine the behavior 
of the eigenfunctions for small values of k. Before tackling the problem at hand we discuss 
a toy model in this appendix. The next appendix discusses the eigenfunctions of the real 
problem. 

The toy model involves the eigenvalue equation 

- dl^Pk + Vip)^k = k^^k (B.71) 

where the potential is a step function 

V{p) = oo,p<0 (B.72) 
= -Vo ,0<p<l (B.73) 
= ,p>l (B.74) 

We want to find the eigenfunctions ipk with continuous spectrum, which clearly starts 
from zero. This potential has features which are similar to the potential Vo(p) in the text. 
The solution for p < 1 is given by 



^fc = ^sin(v/pTT^p) 



< p < 1 (B.75) 



The solution for p > 1 is 



iJk = ^ sm(kp + dik)) (B.76) 



The fc-independent coefficient in (B.76) follows from normalization of the wave function. 
This may be verified by solving the problem with a IR cutoff at some large value of p = L, 
imposing Dirichlet boundary conditions there, and finally taking the limit L — > oo. 
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Even though there is a continuous spectrum starting at = the Green's function for 
the operator on the left hand side of (B.71) exists for generic values of Vq. However for 

VVo = {n+^)7r n = 0,l,2,--- (B.77) 

there is a zero mode solution which vanishes at p = and is regular at p = oo, 

Mp) = ^ sin[(n + l)7rp], p < 1 (B.78) 

'71 2 



1 

P>1 



For these values of the potential, the Green's function diverges. The potential for n = 
above corresponds to a critical point. 

Matching two solutions at p = 1 we get 

A{k) sin(VA;2 + Vo) = sin(A; + e{k)) (B.79) 
A{k)^/k■^ + Vocos{^/k^ + Vq) = kcos{k + e{k)) (B.80) 



leading to 



A{k) = , ^ = - (B.81) 

Oik) = tan-i I ^ 1 - k. (B.82) 



We now examine the behavior of the wavefunctions in the small k limit. For a generic 

2 

potential Vq 7^ ^, the leading order behavior for small k is given by 

A{k) ^ k \ + Oik') (B.83) 

^/Vo COs(VVo) 

0^k) ^ k I^ ^^^J^) -1^+0 (e) (B.84) 

Now consider a function J(p), with finite support. The convolution of J(p) with the 
eigenfunction ip^, in the A; — t- limit becomes 

/»1 /"OO 

Ji, = A{k)J dk sin(VA;2 + 1/op)J(p) + J dk sm{kp + e{k))J{p) (B.85) 

p1 poo 

^A{k)j dksm{^/W^op)J{p) + j dk{kp + e{k))J{p) (B.86) 

~ 0{k) (B.87) 
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Now let us discuss the critical case Vq = In this case we get, upto an overall 
/c-independent constant, 

A{k) ^ I - — + O {k"") (B.88) 
8 

W^-^-^ + 0(A;^) (B.89) 

This leads to 

/I poo 
dk sm{^/k^ + Vop)J{p) + J dk sm{kp + e{k))J{p) (B.90) 

1,2 rl rco i 

dksm{^/¥+Vop)J{p) + dk{l--{kp + k/2f)J{p) (B.91) 
~ Jo + 0{P) (B.92) 

C k dependence of the coefficients : the real problem 

In this appendix we will prove the scaling of quantities like (J7fc — Jo) which appear in 



(4.58). 



The eigenvalue problem we wish to solve is, 

-dlMp) + v{p)Mp) = k'Mp)- (C.93) 

is a continuous parameter which may be chosen to be positive. We will assume some 
simple properties for the potential V{p) near p = oo, i.e., V{p) ~ e'^ as p — oo. Near 
the boundary p = the potential V{p) is assumed to be sufficiently well behaved so that 
the above EOM has a solution of the form, 

V'fc(O) ~ Ap"+ +5p"-. (C.94) 

Where A and B are coefficients of leading and sub leading modes. The potential V{p) 
in the problem of interest (i.e. the equation in a black brane background) satisfies these 
properties. We use the boundary condition, 

A = 0. (C.95) 

The normalization we use is, / dpipk{p)4'k'{p) = ^{k — k'). 



To solve the equation (C.93) for small fc, we divide the range of p in two overlapping 
parts: (i) in the "far" region p ^ 1 we may use V{p) ^ — Vi exp(— p) (ii) in the "near" re- 
gion we can treat k"^ purturbatively. Then we match these two solutions in an overlapping 
region to find a small k solution. 
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In the perturbative "far" region: 

Mp) = m{Mp) + 0{k^)). (C.96) 

Here ipo{p) is the solution of the = equation with the correct boundary conditions at 
p = and A{k) is a constant to be determined. Near p ^ oo, ipo ap + b. The boundary 
condition we choose for the necessary derivative of ipo is such that b = 1. 



For p^ 1 the equation (C.93) may be approximated as, 

dpMp) + (Vi exp(-p) + k^)Mp) = 0. (C.97) 
The real solution of the above equation is given by, 

Mp) = BuJ2rk (2^V^e-p) + 5* J_2.fc {2^V^^ (C.98) 
For p — )■ oo, we have: 

Demanding that the above solution reaches a delta function normalizable sinusoidal solu- 
tion we have, 

BkBl = h:{2ik + l)T{-2ik + l). (C.lOO) 

This condition is needed for the normalization we use and is similar to the choice of 



coefficient in (B.76) in the toy problem considered in the previous appendix. 

In the small k limit we have a matching region (1 ^ po ^ I log(|^)|) where we may 
expand the Bessel functions in small argument and match with the perturbative solution. 
We get at lowest order in fc, 

Bk + Bl = A{k)b (C.lOl) 
2ik{Bk - Bl) = A{k)a (C.102) 

Hence generically for small k, A{k) ~ k ( note 6=1). 

Things are special at the critical point, as a = there. At the critical point Bk = Bl 
and A{k) = 1 ~ o(l). 

Now consider a function J{p) with a compact support such that, 

Jo= [ Jip)Mp)dp, (C.103) 
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is finite. Now generically, 

Jk = J Jip)Mp)dp,^Aik) J J{p)Mp)dp, (C.104) 

~ 0(A;)asA; ^ 0. (C.105) 

In contrast, at the critical point, 

Jk = j Jip)Mp)dp, (C.106) 

^ J Jip)iAik)Mp) + oik')) dp, (C.107) 
~ Jo + o(^^) as A; ^ 0. (C.108) 

The subleading piece comes from the perturbative fc^ corrections. 
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